Electrical Detection of Coherent Nuclear Spin Oscillations in Phosphorus-Doped 

Silicon Using Pulsed ENDOR 
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We demonstrate the electrical detection of pulsed X-band Electron Nuclear Double Resonance 
(ENDOR) in phosphorus-doped silicon at 5 K. A pulse sequence analogous to Davies ENDOR in 
conventional electron spin resonance is used to measure the nuclear spin transition frequencies of the 
31 P nuclear spins, where the 31 P electron spins are detected electrically via spin-dependent transi- 
tions through Si/Si02 interface states, thus not relying on a polarization of the electron spin system. 
In addition, the electrical detection of coherent nuclear spin oscillations is shown, demonstrating 
the feasibility to electrically read out the spin states of possible nuclear spin qubits. 



Control and readout of electron spin states in solids 
are well established down to the single spin level [l|, Q ■ 
In contrast, the readout of nuclear spin states has mostly 
been limited to optical techniques 0, 0]. However, for 
nanostructures not exhibiting luminescence, an electri- 
cal readout scheme is required. Electrical detection of 
nuclear magnetic resonance has been performed on e.g. 
two-dimensional electron gases, identifying the origin of 
the Overhauser held [5[ , and the electrical readout of nu- 
clear spin states has been achieved for this system [f|. 
In phosphorus-doped silicon, McCamey et al. 0] have 
recently demonstrated the electrical readout of nuclear 
spins at 8.6 T. While both these studies employ highly 
polarized spin systems for the readout, we make use of 
a spin-dependent recombination process via Si/Si02 in- 
terface states Q ; this approach does not rely on a polar- 
ization of the electron spin system and thus works un- 
der experimental conditions where the thermal energy 
is much larger than the electron Zeeman splitting [9j. 
We here demonstrate the electrical detection of coher- 
ent nuclear spin oscillations in phosphorus-doped silicon 
using pulsed Electrically Detected Electron Nuclear Dou- 
ble Resonance (EDENDOR) at 0.3 T using X-band fre- 
quencies (10 GHz). Outside the framework of quantum 
computation, pulsed EDENDOR could become a valu- 
able spectroscopic tool for the characterization of point 
defects in semiconductors combi ning the high sensitivity 
of continuous-wave EDENDOR 10] with the reduction 
of the dynamic complexity of the coupled electron and 
nuclear spin systems via pulsed spectroscopy [111 ]. 

The basic principle of pulsed EDENDOR is depicted 
in Fig. [TJ For the 31 P donor in silicon to be investi- 
gated, its electron spin 31 P e with S — 1/2, its nuclear 
spin 31 P n with / = 1/2, and their hyperfine coupling 
give rise to a four-level system. For the electrical readout 
we use a spin-to-charge conversion mechanism based on 
a spin-dependent recombination involving the 31 P donor 
electron and the Pbo center (S — 1/2) at the Si/Si02 in- 
terface [H . Accounting for the two orientations of the Pbo 
spin, we sketch the eight different states for the three in- 



volved spins in panel (i), indicating the occupation of the 
different states by the gray bars. Due to the Pauli prin- 
ciple, spin pairs with antiparallel spins recombine while 
the pairs with parallel orientation are long-lived. There- 
fore, in the steady state only levels associated with par- 
allel orientations of 31 P -Pbo pairs are occupied [9j with 
a probability of 1/4 each. We neglect the polarization of 
the spin system at hrst. 

The presented pulsed EDENDOR experiments are 
based on the Davies pulse sequence for conventional 
pulsed ENDOR [H. At hrst, we describe the EDEN- 
DOR experiments assuming that no recombination oc- 
curs on the timescale of the pulse sequence. The prepa- 
ration microwave (mw) 7r-pulse inverts the populations 
of the levels associated with one of the two 31 P C hyper- 
fine transitions (i). Subsequently, a radiofrequency (rf) 
7r-pulse inverts the populations on one of the 31 P n hyper- 
fine transitions (ii). Since the signal observed in pulsed 
electrically detected magnetic resonance is proportional 
to the fraction [ap] of 31 P e -Pbo pairs with antiparallel 
spins at the end of the pulse sequence [l3j |. a detection 
mw-7r-pulse is applied to electrically readout the nuclear 
spin state (iii) . At the end of this pulse sequence, we ex- 
pect an antiparallel spin fraction of [ap] = 1/2 when an 
rf 7r-pulse has been applied (iv) whereas [ap] = without 
an rf pulse or with an rf pulse far off resonance (iv)*. 
In contrast to conventional ENDOR, where the ENDOR 
intensity is limited by the polarization of the electron 
spin ensemble, the EDENDOR intensity is given by [ap], 
which can be significantly larger. 

For the measurements of EDENDOR, we used a sam- 
ple consisting of a 22 nm-thick Si layer with a nominal 
P concentration of 9 x 10 16 cm~ 3 grown by chemical va- 
por deposition on a 2.5 /im thick, nominally undoped Si 
buffer on a silicon-on-insulator substrate. The epilayer 
is covered with a native oxide. Interdigit Cr/Au con- 
tacts with a periodicity of 20 /im defining an active area 
of 2 x 2.25 mm 2 were evaporated for photoconductivity 
measurements. The sample was mounted with the silicon 
[110] axis parallel to the static magnetic field, illuminated 
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FIG. 1. The pulse sequence of the electrically detected Davies ENDOR consists of a preparation microwave- 7r-pulse, an rf pulse 
of length T p , and a mw detection 7r-pulse after a waiting time t w . This pulse is replaced by a mw detection echo-sequence, 
where the final 7r/2-pulse is phase shifted by 180 degrees with respect to the other pulses to form an effective 7r-pulse. The 
evolution of the spin state populations for an ideal EDENDOR is depicted in the panels in the upper part of the figure. The 
donor electron spin ( 31 P e ), its nuclear spin ( 31 P n ), and the electron spin of the Si/Si02 interface state (Pbo) are symbolized by 
the three arrows in the bracket. The populations of the spin states are indicated by gray boxes. The fraction of antiparallel 
spin pairs at the end of the pulse sequence is given by [ap\. Panels (iii) 1 -(vi) 1 illustrate the evolution of the populations when 
recombination of antiparallel spin pairs during the pulse sequence is taken into account. 



with above-bandgap light, and biased with 100 mV re- 
sulting in a current of «60//A. All experiments were per- 
formed at «5 K in a BRUKER dielectric microwave res- 
onator for pulsed X-band ENDOR. The microwave-pulse 
power was adjusted such that the 7r-pulse length was 
30 ns, corresponding to a microwave B\ -field of 0.6 mT. 
The current transients after the pulse sequence were 
recorded by measuring the voltage drop over a 1.6 kil 
resistor placed in series with the sample. The resulting 
voltage transients were filtered by a cascade of low-pass 
filters (/3dB ~ 4 MHz), amplified, and recorded with a 
digitizer card. The pulse sequence used is shown on the 
top of Fig. [T] To remove the background resulting from 
non- resonant photocurrent transients [Uj , we replace the 
final detection mw-7r-pulse by a detection echo-sequence 
0, [l5| ; by applying a two-step phase cycle [lj| to the last 
microwave 7r/2-pulse (indicated in Fig.[T]by ±x), with the 
phase changed by 180 degrees for every shot, and sub- 
tracting consecutive traces, only resonant changes in the 
photocurrent are detected. To obtain a sufficient signal- 
to-noise ratio, the experiment is repeated with a shot 
repetition time of 800 /is allowing the electron spin sys- 
tem to relax to its steady-state [lg, [l7 1 . The recorded 



current transients are box-car integrated, resulting in a 
charge AQ that is proportional to the antiparallel spin 
pair fraction [ap] at the end of the microwave pulse se- 



quence 



13] 



In a first experiment, we demonstrate that pulsed 



EDENDOR can be used as a spectroscopic method. To 
this end, we keep the length of the rf pulse fixed at 
T p =10/is while sweeping the rf frequency. For the re- 
mainder of this work, the detection echo-sequence with 
equal evolution times of 100 ns is applied after waiting for 
t w =2 /is, c.f. Fig. [TJ Figure [5] (a) shows the integrated 
current transient AQ as a function of the rf frequency 
/ r f for two different magnetic fields, chosen such that 
the microwave pulses are resonant with the high-field 
hyperfine-split 31 P C resonance (upper trace) and with 
the Pbo centers at </=2.0042 (lower trace) [18J. Since 
for non-resonant rf pulses the antiparallel spin content 
is the same for microwave pulses resonant with the 31 P e 
and with the Pbo spins, comparison of the two traces 
allows to identify the 31 P n transitions at frequencies of 
/i=52.25±0.02MHz and / 2 =65.08±0.02MHz. More de- 
tailed frequency scans of the two peaks are shown in 
panels (b) and (c) where the Pbo trace was subtracted 
from the 31 P e trace to remove the non-resonant back- 
ground signal. The measured nuclear transition frequen- 
cies are in good agreement with the expected frequencies 
of 52.34 MHz and 65.19MHz for the hyperfine interaction 
of 117.53 MHz between 31 P C and 31 P n , and the 31 P n nu- 
clear Larmor frequency of 6.076 MHz [lj^. The slightly 
smaller value of the measured frequencies could be at- 
tributed to small deviations of the hyperfine interaction 
at the surface from the bulk value due to local strain 
in the thin- film silicon sample [20| | . To exclude spuri- 
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FIG. 2. (a) Integrated current transient AQ after the Davies 
EDENDOR pulse sequence with T p =10 /is as a function of 
the rf frequency. When the mw pulses are in resonance with 
the 31 P high- field hyperfine line (upper trace) two peaks can 
be observed at / rf =52.25 MHz and / rf =65.08 MHz. These 
peaks do not appear in the lower trace where the mw pulses 
are in resonance with the Pbo spins. More detailed frequency 
scans of the two peaks are shown in panels (b) and (c), where 
the non-resonant background has been subtracted. AQendor 
denotes the peak amplitude. 



ous effects of the intense rf pulse, we repeated the same 
measurement without the preparation mw-7r-pulse. As 
expected, no peaks at the nuclear transition frequencies 
were observed (data not shown). 

To demonstrate the electrical readout of coherent nu- 
clear spin oscillations, we measured the amplitude of the 
EDENDOR signal as a function of the rf pulse length 
T p in a second experiment. To this end, we recorded rf 
frequency sweeps as those in Fig. [5] (b) for different T p 
and fitted the peaks with Lorentzians. In Fig. EH we show 
their amplitudes AQendor a s a function of T p , revealing 
a damped oscillation. We attribute this oscillation to the 
coherent driving of the 31 P nuclear spins. To corroborate 
this interpretation, we measured the oscillation period 
Tpiabi for different rf power levels P r f resulting in a linear 
increase of the Rabi frequency 1 / TRabi with the rf I^-field 
(B 2 ~ ) as shown in the inset of Fig. [3] To compare 
the measured EDENDOR amplitude AQendor with the 
theoretically expected maximum value [op] = 1/2 from 
the considerations depicted in Fig. [TJ we take the ampli- 
tude AQ cc ho of a reference detection echo defined above 
(which corresponds to an effective 7r-pulse) without pre- 
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FIG. 3. Amplitude of the EDENDOR peak AQendor at 
52.25 MHz as a function of the rf pulse length T p . The os- 
cillation shows the coherent driving of nuclear spin motion 
which is damped on the timescale of the recombination of the 
31 P e -Pbo spin pairs. The Rabi frequency scales linearly with 
the square root of the rf power, as shown in the inset. 



ceding mw and rf pulses under the same experimental 
conditions as a measure for the maximum experimen- 
tally achievable value of [ap\. We define the EDENDOR 
contrast as AQENDOR/AQccho, resulting in a contrast of 
0.07 at T p =10^s. 

The evolution of the spin system during the Davies 
EDENDOR sequence deviates in several aspects from the 
ideal situation discussed above. These are primarily the 
recombination that occurs during the rf pulse and the 
polarization and relaxation of the electron and nuclear 
spins. A full simulation of the experimentally observed 
data will require incorporation of these dynamic effects 
e.g. in a system of combined rate equations [ljjj]. To gain 
a first physical picture, we here separately discuss the ef- 
fects of the recombination time r ap characteristic for an- 
tiparallel electron spin pairs, of the recombination time 
T p characteristic for the parallel spin pairs, and of the po- 
larization of the nuclear spin system. To independently 
quantify the recombination time r ap , we performed an 
inversion-recovery type of experiment 17]. The observed 
decay can be described with a stretched exponential with 



a time constant of 1 [is, which we identify with r ap , and 
an exponent of ~ 0.5. (The recombination time r p is not 
readily accessible; we assume that it is much longer than 
T ap due to the weak spin-orbit coupling in Si and that 
is shorter than the shot repetition time used.) Thus, af- 
ter an rf pulse with a length of T p =10/is, we expect an 
EDENDOR amplitude of 0.3 • AQ ocho . In addition, the 
finite excitation bandwidth of the rf pulse only excites 
about a fraction of 0.6 of the nuclear spin resonance line 
with a FHWM of 0.12 MHz. This reduces the expected 
EDENDOR contrast to ssO.18, which agrees within bet- 
ter than a factor of 3 with the observed EDENDOR con- 
trast despite the crude approximations involved. 
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As Fig. [3] demonstrates, even for times T p much longer 
than T a p we observe an EDENDOR oscillation amplitude. 
This can be qualitatively explained in terms of the evo- 
lution of the spin system depicted in Fig. [T] To esti- 
mate the effect of longer rf pulses T p , we hypothetically 
divide T p into an ideal ir pulse without recombination 
during the pulse, a time period t, during which the an- 
tiparallel spin pairs recombine, and an additional ideal rf 
pulse. Starting from the spin system after the first ideal 
7r rf pulse (iii) , we consider the situation after a time t 
(r p > t> r ap ) when all antiparallel spin pairs have recom- 
bined leaving unoccupied donor states. These are refilled 
with electrons that form new parallel spin pair states with 
equal probability for each of the four parallel spin-pair 
states (iv) 1 , since the formation of new spin pairs is not 
spin-dependent Q . The additional ideal rf pulse induces 
nuclear spin oscillations that can be measured by a detec- 
tion echo because of the remaining population difference. 
For example, after an additional rf 7r-pulse we expect 
[ap] = 1/4 as shown in (iv) 1 -(vi) 1 , whereas an additional 
2ir rf pulse would result in [ap] = 1/2. Therefore, the 
expected amplitude of the oscillation is reduced by a fac- 
tor of 2 compared to the oscillations expected according 
to (ii)-(iv). Motivated by the above considerations, we 
tentatively attribute the experimentally observed persis- 
tence of the oscillation and its amplitude at pulse lengths 
T p S> T ap to the substantial difference in the involved re- 
combination times. These oscillations are expected to 
decay on the timescale of t p , which would be accessible 
by EDENDOR according to this line of argument. 

In analogy to conventional Davies ENDOR 12], the 
ideal EDENDOR pulse sequence transfers the electron 
spin polarization to the nuclei. This transfer of polariza- 
tion allows to readout the nuclear spin state after wait- 
ing times t w much larger than the recombination times, 
limited only by Ti n . Indeed, we observe an EDENDOR 
signal for t w > 4 ms with an EDENDOR contrast of 0.01 
using a shot-repetition time of 15 ms. This would allow 
to measure Ti n by recording the EDENDOR contrast as 
a function of t w . 

In conclusion, we have demonstrated the electrical de- 
tection of coherent nuclear spin oscillations with pulsed 
EDENDOR employing Si:P as a model system. We use 
a spin-dependent recombination process via Si/Si02 in- 
terface states enabling the nuclear spin readout without 
making use of a polarization of the electron spin sys- 
tem. During the readout process, the electron is removed 
from the donor, which drastically changes the Hamilto- 
nian describing the spin system. For applications of this 
approach as quantum memory 0, Hi} , such effects might 
be suitably suppressed by decoupling the 31 P e -Pbo read- 
out pair using local gate electrodes, where an applica- 
tion of an electric potential keeps the donor in its para- 
magnetic state for a long time. On the other hand, it 
would be interesting to study the effects of removing and 
repopulating the electron on the nuclear spin state in 



detail, e.g. with respect to the relaxation and decoher- 
ence of the nuclear spin system. Furthermore, pulsed 
EDEDNOR gives access to the recombination time of 
parallel spin pairs that has hitherto eluded unambiguous 
experimental determination as well as the nuclear relax- 
ation time 0. The high sensitivity of the electrical de- 
tection of magnetic resonance and the widespread avail- 
ability of X-Band spectrometers makes the presented ap- 
proach to pulsed EDENDOR an attractive spectroscopic 
method for the characterization of defects in semicon- 
ductor nanostructures. Hereby, hyperfine and superhy- 
perfine interactions can be investigated, which allow to 
obtain quantitative information about the electron wave- 
function and the specific environment around the defect. 
Since it does not rely on polarization of the spin sys- 
tem, this approach could also be particularly attractive 
for room temperature EDENDOR spectroscopy, e.g. of 
endohedral fullerenes such as N@C60 [22J. 

The work was supported by DFG (Grant No. SFB 631, 
C3). 
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